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Abstract. We consider the relativistic electron-positron field interacting with 
itself via the Coulomb potential defined with the physically motivated, positive, 
density-density quartic interaction. The more usual normal-ordered Hamilton- 
ian differs from the bare Hamiltonian by a quadratic term and, by choosing 
the normal ordering in a suitable, self-consistent manner, the quadratic term 
can be seen to be equivalent to a renormalization of the Dirac operator. For- 
mally, this amounts to a Bogolubov-Valatin transformation, but in reality it is 
non-pcrturbative, for it leads to an inequivalent, fine-structure dependent rep- 
resentation of the canonical anticommutation relations. This non-perturbative 
redefinition of the electron/positron states can be interpreted as a mass, wave- 
function and charge renormalization, among other possibilities, but the main 
point is that a non-perturbative definition of normal ordering might be a useful 
starting point for developing a consistent quantum electrodynamics. 



1. Introduction 

In relativistic quantum electrodynamics (QED) the quantized electron-positron 
field ty(x), which is an operator- valued spinor, is written formally as 

(1) 9(x) := a(x) + b*(x), 

where a(x) annihilates an electron at x and b*(x) creates a positron at x. (We 
use the notation that x denotes a space-spin point, namely x — (x, a) G T — 
R 3 x {1,2,3,4} and J d 3 x denotes integration over R 3 and a summation over the 
spin index.) More precisely, we take the Hilbert space, S) of L 2 (R 3 ) spinors, i.e., 
f) = L 2 (R 3 ) <g) C 4 . To specify the one-electron space we choose a subspace, of 
and denote the orthogonal projector onto this subspace by P+. The one-positron 
space is the (anti-unitary) charge conjugate, Ci}_, of the orthogonal complement, 
$j- with projector P_. In position space, (C(ip))(x) := i/Ja^^x), whereas in 

momentum space (C(ip))(p) := if3at2'ip(—p). For later purposes we note explicitly 
that the Hilbert space 9) can be written as the orthogonal sum 

(2) £=io+ei}_. 

If /„ is an orthonormal basis for S) + with v > and an orthonormal basis for 
with v < then a(x) := Y%Lo a (fv)fv( x )> and b(x) := J2u=-oo b (f")f>^( x )^ where 
a*(f) creates an electron in the state P+/, etc. (For further details about the 
notation see Thaller Q , and also Helffer and Siedentop Q| and Bach et al ^, llj. 

1 © 1999 by the authors. This paper may be reproduced, in its entirety, for non-commercial 
purposes. 
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For free electrons and positrons, Sj+ and io_ are, respectively, the positive and 
negative energy solutions of the free Dirac operator 



in which a and (3 denote the four 4x4 Dirac matrices and p = — iV. The number mo 
is the bare mass of the electron/positron. Perturbation theory is defined in terms 
of this splitting of Sj into $)+ and $)-. The electron/positron lines in Feynman 
diagrams are the resolvents of Do split in this way. As we shall see, this splitting 
may not be the best choice, ultimately, and another choice (which will require a fine- 
structure dependent, inequivalent representation of the canonical anticommutation 
relations) might be more useful. 

The Hamiltonian for the free electron-positron field is 



with this choice of The symbol : : denotes normal ordering, i.e., anti-commuting 
all a*,b* operators to the left (but ignoring the anti-commutators). 

In this paper we investigate the effect of the Coulomb interaction among the 
particles, which is a quartic form in the a # ,6 # operators. (# denotes either a 
star or no star.) The normal ordering will give rise to extra constants and certain 
quadratic terms. Our main point will be that by making an appropriate choice of fi + 
that is different from the usual one mentioned above, we can absorb the additional 
quadratic terms into a mass, wave function, and charge renormalization. This 
choice of fj + has to be made self consistently and, as we show, can be successfully 
carried out if some combination of the fine structure constant a (=1/137 in nature) 
and the ultraviolet cutoff A is small enough (but 'small' is actually 'huge' since the 
condition is, roughly, a log A < 1). 

The physical significance of our construction is not immediate, but it does show 
that certain annoying terms, when treated non-perturbatively, can be incorporated 
into the renormalization program. It is not at all clear that what we have called 
wave function renormalization, for example, really corresponds to a correct inter- 
pretation of that renormalization, or even what the true meaning of wave function 
renormalization is. It could also be interpreted as a renormalization of Planck's 
constant. Likewise, the charge renormalization we find is not mandatory. But our 
main point is that the effect is mathematically real and is best dealt with by a 
change of the meaning of the one-particle electron/positron states. It should be 
taken into account in the non-perturbative QED that is yet to be born. Of course, 
there are other renormalization effects in QED that we do not consider. In par- 
ticular, the magnetic field has not been included and so we do not have a Ward 
identity to help us. Our result here is only a small part of the bigger picture of a 
proper, non-perturbative QED. 

The last section of this paper contains a brief discussion of some possible inter- 
pretations of our findings and the interested reader is urged to look at that section. 

Our starting point is the unrenormalized ('bare') Hamiltonian for a quantized 
electron/positron field interacting with itself via Coulomb's law, namely 



(3) 



D = a ■ p + m Q (3, 




(4) HP 



bare 



d 3 x : V*(x)D if>(x) : +- / d 3 x / d 3 yW(x,y) : **(a;)tf (a) :: 
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where W is a symmetric (W(x, y) = W(y,x)) interaction potential. The case of 
interest here is the Coulomb potential W(x,y) — £ CT , T |x — y| _1 and regularized 
versions thereof. We will refer to the first term on the right side of (||) as the 
kinetic energy T and to the second term as the interaction energy W bare . 

Our choice of the product of the two normal ordered factors in the interaction, 
namely : p{x) : : p{y) :, is taken from the book of Bjorken and Drcll jj| (Eq. (15.28) 
without magnetic field). It is quite natural, we agree, to start with the closest 
possible analog to the classical energy of a field. Of course, this term is only partly 
normal ordered in the usual sense and therefore it does not have zero expectation 
value in the unperturbed vacuum or one-particle states. It is positive, however, as 
a Coulomb potential ought to be. 

The fully normal ordered interaction (which is not positive) is defined to be 

(5) W r Q cn = ^Jd 3 xJ d 3 yW(x,y) : V* (x)V* (y)V (y)V (x) : . 

Bear in mind that this definition entails a determination of the splitting (0) of 9). 
We introduce a normal ordered ('dressed') Hamiltonian 

(6) H^ n := J d 3 x : ^*{x)D z , m ^{x) : +W rcn °. 

(Again, this depends on the splitting of Sj.) The operator Dz, m differs from Do m 
two respects: 

(7) D Ztn = Z- 1 (a-p + mp), 

where Z < 1 and m > tuq. We call mo the bare mass and m the 'physical' mass — 
up to further renormalizations, which we do not address in this paper. The factor Z 
is called the 'wave function renormalization' because the 'wave function renormal- 
ization' in the standard theory is defined by the condition that the electron-positron 
propagator equals Z times a non-interacting propagator, which is the inverse of a 
Dirac operator. To the extent that W™ n can be neglected, the propagator with just 
the first term on the right side of @) would, indeed, be Z times a Dirac operator 
with mass m. One could also interpret Z in other ways, e.g., as a renormalization 
of the speed of light, but such a choice would be considerably more radical. Some 
other interpretations are discussed in the last section. 

Our goal is to show that by a suitable choice of Z, m, and the electron space 
S)-\-, we have - up to physically unimportant infinite additive constants - asymptotic 
equality between H barc and H^ cn for momenta that are much smaller than the cutoff 
A. 

(8) H baro = H™ n (for small momenta) 

The additional normal ordering in (^|) will introduce some quadratic terms and 
it is these terms to which we turn our attention and which we would like to identify 
as renormalization terms. Whether our picture is really significant physically, or 
whether it is just a mathematical convenience remains to be seen. What it does do 
is indicate a need for starting with a non-conventional, non-perturbative choice of 
the free particle states, namely the choice of 
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2. Calculation of the New Quadratic Terms 
The canonical anti-commutation relations are 

{«(/), a(g)} = {a*(f),a*(g)} = {«(/), b(g)} = {a*(f),b*(g)} 
( ' ={a*(f),b(9)} = Hf),b*(g)} = 0, 

and 

(10) {a(f), a*(g)} = (/, P + g), {&*(/), 6(g)} - (/, P_ ff ). 

Formally, this is equivalent to 

{a(x),a(y)} = {a* (x) , a* (y)} = {a(x),b(y)} = {a* (x) , b* (y)} 
( } ={a*(x),b(y)} = {a(x),b*(y)} = 0, 

and 

(12) {a(x),a*(y)} = P+(x,y), {b*(x),b(y)} = P-(x,y), 

where P±(x,y) is the integral kernel ol the projector P±. 

In order to compare the bare Hamiltonian H barc with the renormalized one H^ on 
we must face a small problem; the momentum cutoff A, which is necessary in 
order to get finite results, will spoil the theory at momenta comparable to this 
cutoff. Thus, our identification of the difference as a wave function and mass 
rcnormalization will be exact only for momenta that are small compared to A. We 
will also require (although it is not clear if this is truly necessary or whether it is 
an artifact of our method) that aln(A/m ) is not too large. 

The potential energy W bare clearly has positive singularities in it if the unrcg- 
ularized Coulomb interaction W(x,y) = |x — yl -1 ^,- is taken. We therefore cut 
off the field operators for high momenta above A. (Another choice would be to cut 
off the Fourier transform of |x — y| _1 at A and the result would be essentially the 
same; our choice is motivated partly by computational convenience.) Our cutoff 
procedure is to require that 

(13) *(ar) := (2tt)- 3/2 / §(p,cr)e- J P- x dp 

J{ P m 3 | | P |<A} 

holds, with A > 0. We could also take a smooth cutoff without any significant 
change. This regularization keeps the positivity of W bare . To make all terms finite 
we also need a volume cutoff for one of the difference terms. However, the volume 
singularity only occurs as an additive constant, which we drop. The renormalized 
Hamiltonian is free of infrared divergent energies and so does not need a volume 
cutoff. 

Next, we calculate the difference 
(14) 

K := W baro -Wr = f J d 3 x J d 3 yW(^y)[a*(x)a(y)P + (x,y)+a*(x)b*(y)P + (x,y) 

+ b(x)a(y)P + (x, y) - b* (y)b(x)P + (x, y) 
- a*(y)a(x)P- (x, y) + P_ (x, y)P+ {x, y) + a(x)b(y)P- (x, y) 

+ b*(x)a*(y)P-(x,y) + b(x)*b(y)P.(x, y)] . 
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Thus, (Q) can be rewritten as 

(15) H barc = [ d 3 x : 9*(x)D ^(x) : 



The first two terms on the right side are one-particle operators (i.e., quadratic 
in the field operators). Let us call their sum A, which we write (formally, since A 
is a differential operator) as 

(16) A= / d 3 x J d 3 yA{x,y) : V*(x)V{y) : . 

The last term is a cutoff dependent constant, which happens to be infinite. (One 
can show that this term is positive and, if we put in a momentum cutoff A and 
volume (infrared) cutoff V, this integral is proportional to VA 4 .) 

We can make A positive by choosing the normal ordering appropriately. That is, 
we choose P+ to be the projector onto the positive spectral subspace of the operator 
A. Thus, we are led to the nonlinear equation for the unknown A 

(17) A = D + -R, 
where the operator R has the integral kernel 

,,o\ t>( \ P +( X >V) - p -{x,y) (sgnA)(x,y) 
18 R(x,y):= : : = — : : — , 

l x -y| l x -y| 

with sgni being the sign of t. 

For physical reasons and to simplify matters we restrict the search for a solu- 
tion to ( |l7| ) to translationally invariant operators, i.e., 4x4 matrix- valued Fourier 
multipliers. Moreover, we make the ansatz 

(19) A(p):=o.«p5i(|p|)+/350(|p|) 

with real functions g\ and go, where lu p is the unit vector in the p direction. In 
other words, we try to make A look as much as possible like a Dirac operator. With 
this ansatz we have sgnA(p) = A(p)/(<7i(|p|) 2 + ^o(IpI) 2 ) 1 ^ 2 - Thus, recalling (|l3|), 
(0) can be fulfilled, if 

(20) 0O(|P|) :=m + ^2 / d % ~~i2 f 7\ h/r~n h2U/2 

4tt 2 J| q |<A |p - q| 2 (gi(|q|) 2 + go(|q|) 2 ) 1/2 

r<r\\ n \\ II, ^ f w P • w q ffi(|q|) 

(2D M) := |P| + ^ j / q|<A ^|p_ q | 2 ( 5l( | q |) 2+ff0 (| q |) 2 ) 1 / 2 

is solvable. Note that the bare mass mo appears in (pTj). 

3. Determination of the Dressed Electron 

We will find a solution of the system (pTj|), ( pl| ) by a fixed point argument. To 
this end we first integrate out the angle on the right hand side. Setting u := |p| 
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and v :— |q| we get 

(22) g (u) = mo 

(23) 9i(u)=' 
where, for z > 1, 
(24) 
and 
(25) 



a 
27 



dv- 



/ 1 /u v 



9o{v) 



a 
2^ 



dv- 



[gi{v? + goiv) 2 ) 1 ' 2 



n ( \ 1 i z + 1 
Q (*) = -log — - 



M ,) = |lo g i±l-l. 



The case of zero bare mass is particularly easy. We can choose go = 0, in which 
case 51 is obtained by integration. We get 



9i W 



— I dv-Qt (-(- + - 



(26) 



2tt 

a 



v 



v 



=u + — — u / dvv 



2tt 

a 

— i 

2tt 



A/u 



A 

6 it 



log 



2 
3 

(A/« 



log|(A/ u ) 2 
1 



i((l/« + «)/2) 
1 



1 



i(A/«) 



(3 + (A/ u ) 2 



(A/u) - 1 

The function behaves asymptotically for small u or large A as gi(u) = 
^:|ulog(A/u). Unfortunately, because of the logu, the operator A can never 
look like the renormalized Dirac operator in (Q) . This asymptotic expansion can be 
seen either from (^6|) or else by noting that the integrand is a continuous function, 
except for v = 1, that decays at infinity as 4/(3u). This follows from the large v 
expansion 



h((l/v + v)/2)= -v-^—v 



+ 0(v~ 5 ). 



(27) 

For positive bare masses we solve the system (|22j) and (|23|) by a fixed point 
argument. To this end we define the following set of pairs of functions, for e, S > 0. 

(28) S e ,s :={g= (go, Si) |V ue[0 ,A]g(u) G [m ,(l + 8)m ] x [tt, (1 + e)u]}. 

Note that with the metric generated by the sup norm S e g is a complete metric 
space. Next define T : S e ^s —* S e ^s by the right hand side of (p3|) and (|22), i.e., 



29 T g u =m + - / cfo-Qo - ~ + ~ 
(30) r 1 (g)(i»)=ti+^ : / A d^Q 1 fI(^ + ^ 



9o(v) 



(9i(v) 2 + goiv) 2 ) 1 / 2 
bgi(v)*+go{v)*y/*' 



Lemma 1. if Y } defined by 

(3i) y 



a arsinh(A/TOo)/7r, 



safe/?es y < 9/50, if e > 50Y/(9 - 50Y), and if 5 > Y/(l - Y), then T maps S e , 5 
into S e s ■ 
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Note that arsinh(x) = log(x + \Jx 2 + 1), i.e., Y grows logarithmically in A/mo- 



Proof. Obviously, To(g)(u) > too and Ti(g)(u) > u. 

To bound go we return to ( po|) and note the bound (which holds in S e> $) 

'| q |<A ip-4r(|q| 2 



ffo(IPl) < mo + 7^7 / dq 1 — — r 
4tt 2 7i q i<A |p - q|" 



(32) 



<mo + 



4tt 2 



rfq 



1 



(1 + 6)m 



|q|<A 



q 2 (|q| : 



TO o[l H — (1 + 6) arsinh(A/mo)]- 



The second inequality holds because the convolution of two symmetric decreasing 
functions is symmetric decreasing. 

Next we turn to gx, which is a bit more complicated. We split the integration 
region in @ into A = {q| |q| < 2|p| and |q| < A} and B = {q| |q| > 2|p| and |q| < 
A}. In region A we use w p -w q < 1 and hence the contribution to g\ from this region 
is bounded above by 



(33) 



a 
4^2 



dq- 



1 



(l + e)2|p| 



< 2|p|-(l + e) arsinh(A/m )], 



'|q|<A |p-q| 2 (|q| a +"^) 1/a 

for the same reason as in ([32]). 

In region B we use that (|p| 2 — |q| 2 ) 2 > 9|q| 4 /16. We also note, for the integration 
over B, that we can take the mean of the integrand for q and — q. In other words, 
we can bound this term as follows: 



(34) 



<- 



a 
4^2 

Q 



9i(|q|) 



b q |p-q| 2 (.9i(|q|) 2 +5o(|q|) 2 ) 1 / 2 



rfq 



|p-q|- 



P + q- 



(5i(|q|) 2 +.9o(|q|) 2 ) 1 / 2 



<SL [ rla 4|p||q| ( 1 + e )lql 

-8W B q (|p| 2 -|q| 2 ) 2 (|q| 2 



2)1/2 



32a(l + e) 



, J 1_ 

B q |q| 2 (|q| 2 +r, 



.2)1/2 



< 



9tt 



p| arsinh(A/mo) 



Thus, we obtain the bound 



(35) 



50a 



ffi(IPl) < |P|{1 + -^-(l + e) arsinh(A/m )}. 

97T 



□ 



Theorem 1. J/Y := aarsinh(A/TO )/7r < 9/50, if (2 + e + S)Y < 1 and i/ e,<5 
satisfy the conditions of Lemma i/ien T is a contraction. 

Proof. Thanks to Lemma |l| we only need to establish the contraction property. We 
first note that for positive real numbers x, y, x, y 



(36) 



(a;2 + y2)l/2 (J2 + £2)1/2 



< 



£2 _|_ ^2 



((x - x)' + (y - yY) 



- f^2U/2 
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where (£,77) is some point on the line between (x,y) and (x, y). 
Thus we get 

\[T (g)-T (~g)}(u)\ + \[T 1 (g)-T 1 (g)}(u)\ 



-2tt 



dv — 



1 fu 



to o \ - -l- - 

2 \j) it 



V2 W + u 



<■ 



(37) 



: 2vr 



< — 
-2vr 



dv — 



1 /u 



x [(flo(f) -So(«)) 

/ dw- Qo ( r - + - 

Jo u L \2\u it 

\2 Vw u/ 



9o(v) 



9o( v ) 



V9o(v) 2 +gi(v) 2 ^goiv) 2 +g 1 ( v y 
9i(v) 9i(v) 



V9o(v) 2 +gi(v) 2 ^goiv) 2 ThJvj 



2 \v ' uV/ v 2 + m 2 . 



+ (gi(v) -gi{v)) 

(l + e)v 



21 1/2 



(l + <J)m 



lis — el 



< 



"2tt 



dv-Qo 
u V 2 



1 /it v\\ 
2\v + u)) 



(1 + e)t> + (1 + 5)m 



<(2 + e + S)y||g-g|| 

where Y = a arsinh(A/mo)/7r as in Lemma [p. In the last line we have simply noted 
that the integrals are smaller than the corresponding integral in (^). □ 

Corollary 1. The map T has a unique fixed point, ifY< 1/7. 

Proof. Take e = 50F/(9 - 5QY) and 8 = Y/(l - Y) so that T maps S e ,s into itself. 
Then, if Y < Yq, the contraction condition will be satisfied. □ 

For a = 1/137, i.e., the physical value of the fine structure constant - the 
condition Y < 1/7 is fulfilled for A/771 < e 18 . 

4. Properties of the Renormalized Hamiltonian 

What we have shown up to now is that the bare Hamiltonian in (^) is equivalent, 
apart from some additive constants, to a renormalized Hamiltonian. This renormal- 
ized Hamiltonian has the form (Q), but the quadratic term is only approximately 
the one given in (Q). What takes the place of Dz, m is given in © and ©. We 
see immediately from (^2|), (|2^) that the solution is a continuous pair of functions. 
We also see that as long as mo is not zero the functions go and g\ are positive and 
behave properly for small |p|, i.e., go is constant and g\ is proportional to |p|. To 
relate this to Dz t m we first factor out the constant lim u ^ +(5i( u ) l u ) an d call this 
Z^ 1 . As is evident from (|2l|), this constant 1/Z is larger than one — as it should 
be. 

The next thing to verify is that Zgo(0) is bigger than mo, since this is the 
renormalized mass m that appears in (|t]). In other words, we have to verify that 



(38) 



m 



7770 



^)>lim^M 

777 P 



z-\ 
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We shall, in fact, prove more than this: 



Theorem 2. Assuming that Y < 1/7, the unique solution to the equation ( P0|) and 
mentioned has the property 



go(lpl) > gjM) 

m |p| 



(39) 
for all p ^ 0. 

Proof. Define S,5 je to be the subset of Ss, e on which 
(40) 



go(lpl) > gjM) 



m |p| 

for all p ^ 0. If we can show that T also leaves Sg,e invariant, then we have shown 
the wanted inequality ( |39"|) on the solution in Sg >e (by the uniqueness of the solution 
on Ss t e) except for the possibility that equality also can occur. This is so, since we 
can apply the fixed point argument not only to S$ ie but also to Ss, e - 

Showing that (fit)]) holds is - according to ( p2[ ) and ( p3|) - equivalent to showing 
that 



(41) 



v 

dv-Q 

u 



1 fu 

2 \v 



g (v)/m 



(9i(v) 2 

A 



goiv) 2 ) 1 ' 2 



> 



dv-Qi ( \ 
u \ 2 



9i(v)/u 



holds. Now, using the Inequality (40) and the fact that the factors with the roots 
are monotone functions go and g\ respectively, it is enough to show that 



(42) 



dv-Qo (-(- + - 



1 



(ml+v 2 ) 1 / 2 

-A 



> 



dv- 



1 /u v\ 

2 \v + u) 



v/u 



{ml 



v 2y/2 



holds. To proceed, we now compare the integrands pointwise. Using the explicit 
expressions ( pi| ) and ( p5| ) for Qo and Q±, and with t = v/u, we will have shown 
that the integrand of the left hand side pointwise majorizes the one on the right 
hand side of @ if 

2t 



(43) 



log 



t+l 



t- 1 



< 



l* 2 -l| 



holds for all t =/= 1. By symmetry, we only have to consider (^3|) for t > 1. To this 
end we exponentiate ( [43| ) 

t+l ( 2t 

— <eXp {#=1, 
and expand the exponential up to second order. (Note that this gives a lower bound 
on the exponential because the argument of the exponential function is positive.) 
I.e., it suffices to show that 



(44) 



t + l 



< 1 



2t 



t- 



t-1 t 2 
which follows by direct computation. 



(t 2 
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Having established (||) we see that ( fi4| ) indeed gives strict inequality in (|J) for 
the unique fixed point in Sg e . □ 



5. Interpretations of our Results 

What we have done is start with the 'bare' Hamiltonian H barc in (Q), in which 
the interaction is the closest analog to the classical electrostatic energy of a field. 
After some analysis, we found in section 3 that Z and m could be uniquely chosen 
so that for momenta much less than the ultraviolet cutoff A, H barc = H^ cn plus a 
well defined infinite constant. Of course the field W in the two cases is different. 
Formally, they differ by a Bogolubov transformation, but in fact an inequivalent, 
a dependent, representation of the CAR is needed. In section 4 we showed that 
not only is Z~ x m > mo but more is true, namely m > m , and this is comforting 
physically. 

We emphasize that we have not 'integrated out' any field variables. Our new 
Hamiltonian HJj en is the same as the original one — on a purely formal level. It 
is suggestive, nevertheless, that a good deal of the electrostatic energy has now 
been incorporated into the leading term J d 3 x : ^* (x)Dz, m ^(x) : and that the 
remaining interaction is somehow less important than the original one. It is, after 
all, normal ordered, which means it vanishes on one-electron states, if we define 
such states by ^/*|0), where |0) is the vacuum of the new ^. While this makes 
sense perturbatively, it is, however, misleading because the new vacuum (the state 
of lowest energy) is surely not the obvious choice |0). 

If we drop the new interaction term W™ n we are left with J d 3 x : ^*(i)-Dz, m $(a;) : 
as our Hamiltonian. Unfortunately this is not the Hamiltonian of a Dirac operator 
(even at low momenta) because of the factor Z~ x . We have called this a wave 
function renormalization, but that is not really in the spirit of renormalizing the 
one-electron states (which is what is usually done) and is, instead, a renormal- 
ization of an operator. One school, (Kallen ||) thinks it is proper to speak of a 
renormalized by incorporating a factor Z -1 / 2 into but this changes the anti- 
commutation relations! Note that this formulation requires renormalizing the bare 
mass mo to go(0)/Z and renormalizing the bare charge e to e/Z. 

Another point of view is to regard Z _1 m as the physical mass and to change 
Z p into p by a unitary transformation (which is nothing other than a change of 
length scale). This has the disadvantage of changing the speed of light or Planck's 
constant. It also would mean a different scale change for particles of different mass, 
e.g., muons. As opposed to Kallen's procedure, this requires renormalizing the mass 
from mo to go(0) only, but there is no need for charge renormalization. 

Another possibility is to bring out the factor Z^ 1 as a multiplier of the whole 
Hamiltonian, which would mean changing the fine structure constant to aZ, i.e., 
a charge renormalization but now from e to e/Z 1 / 2 only. The obvious problem 
here is that since the Hamiltonian is the generator of time translation, this means 
a change of the time scale (which, again, depends on the particle in question). 

Doubtless, different people will have different opinions about these matters. We 
do not wish to commit ourselves to any point of view. But it is our opinion that 
the construction of W£ n is a significant piece of the puzzle of constructing a non- 
perturbative QED. 
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